ValueFolderTemp='/sps/isc/rneveu/EtudeIRM/Nature/ResultsJournal/';
ValueTime=2;
ValueFilter=1;
n1=32;
ValueFile=['FileAnalysis' num2str(ValueTime) 'Filtered' num2str(ValueFilter) 'Journal' num2str(n1) '.mat'];
load([ValueFolderTemp ValueFile]);
n2=1;
TestTemp=0;
while n2<=size(VecResultsArticle2,1)&&TestTemp==0
    if VecResultsArticle2(n2,4)>0
        TestTemp=1;
    end
    n2=n2+1;
end
VecValueTemp=VecResultsArticle2(n2-1:end,:);
VecValueTemp2=[VecValueTemp(:,4) VecValueTemp(:,3)./VecValueTemp(:,4)*100];
VecValueTemp3=VecValueTemp2(size(VecValueTemp2,1)-30:size(VecValueTemp2,1),:);
VecValueTemp4=VecValueTemp2(1:size(VecValueTemp2,1)-31,:);
VecValueTemp3=[(VecValueTemp3-repmat(mean(VecValueTemp3),size(VecValueTemp3,1),1)).^2 abs(VecValueTemp3-repmat(mean(VecValueTemp3),size(VecValueTemp3,1),1))];
VecValueTemp4=[(VecValueTemp4-repmat(mean(VecValueTemp4),size(VecValueTemp4,1),1)).^2 abs(VecValueTemp4-repmat(mean(VecValueTemp4),size(VecValueTemp4,1),1))];
VecValueTemp5=1:size(VecValueTemp3,1);
VecValueTemp6=1:size(VecValueTemp4,1);
% Comparison of the periods 1945-1989 and 1990-2020 for the articles and the ratio of original articles published in Nature
VecResults=[];
for n1=1:size(VecValueTemp3,2)
    [ValueResults,ValueResults2]=ttest2(VecValueTemp3(:,n1),VecValueTemp4(:,n1));
    [ValueResults3,ValueResults4]=ranksum(VecValueTemp3(:,n1),VecValueTemp4(:,n1));
    [ValueResults5,table,stats]=anovan(VecValueTemp3(:,n1),{VecValueTemp5'},'continuous',1,'model','full','display','off');
    [ValueResults6,table,stats2]=anovan(VecValueTemp4(:,n1),{VecValueTemp6'},'continuous',1,'model','full','display','off');
    VecResults=[VecResults;mean(VecValueTemp3(:,n1)) sqrt(var(VecValueTemp3(:,n1))) mean(VecValueTemp4(:,n1)) sqrt(var(VecValueTemp4(:,n1))) ValueResults2 ValueResults3 NaN ValueResults5 stats.coeffs' NaN ValueResults6 stats2.coeffs'];
end

                    